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Abstract 

We present a numerical simulation of the double slit interference experiment realized by F. 
Shimizu, K. Shimizu and H. Takuma with ultracold atoms. We show how the Feynman path 
integral method enables the calculation of the time-dependent wave function. Because the evo- 
lution of the probability density of the wave packet just after it exits the slits raises the issue of 
the interpreting the wave/particle dualism, we also simulate trajectories in the de Broglie-Bohm 
interpretation. 
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I. INTRODUCTION 



In 1802, Thomas Young (1773-1829), after observing fringes inside the shadow of playing 
cards illuminated by the sun, proposed his well-known experiment that clearly shows the 
wave nature of light. ^ He used his new wave theory to explain the colours of thin films 
(such as soap bubbles), and, relating colour to wavelength, he calculated the approximate 
wavelengths of the seven colours recognized by Newton. Young's double slit experiment is 
frequently discussed in textbooks on quantum mechanics.'^ 

Two-slit interference experiments have since been realized with massive objects, such 
as electronsj ^ * ^ * ^ * ^ * neutrons,'^ cold neutrons,'^ atoms,'^ and more recently, with coherent 
ensembles of ultra-cold atoms,'^^'^ and even with mesoscopic single quantum objects such 
as Ceo and C70.'™' 

This paper discusses a numerical simulation of an experiment with ultracold atoms re- 
alized in 1992 by F. Shimizu, K. Shimizu, and H. Takuma.'^ The first step of this atomic 
interference experiment consisted in immobilizing and cooling a set of neon atoms, mass 
m = 3.349 X 10~^^kg, inside a magneto-optic trap. This trap confines a set of atoms in a 
specific quantum state in a space of ~ 1 mm, using cooling lasers and a non-homogeneous 
magnetic field. The initial velocity of the neon atoms, determined by the temperature of 
the magneto-optic trap (approximately T = 2.5 mK) obeys to a Gaussian law with an 
average value equal to zero and a standard deviation = \J~^^ — 1 m/ s; ks is Boltzman's 
constant. 

To free some atoms from the trap, they were excited with another laser with a waist of 
30 /im. Then, an atomic source whose diameter is about 3 x 10"^ m in and 10^^ m in the 
z direction was extracted from the magneto-optic trap. A subset of these free neon atoms 
start to fall, pass through a double slit placed at li = 76 mm below the trap, and strike 
a detection plate at I2 = 113 mm. Each slit is 6 = 2 fim wide, and the distance between 
slits, center to center, is d = 6 /im. In what follows, we will call "before the slits" the space 
between the source and the slits, and "after the slits" the space on the other side of the 
slits. The sum of the atomic impacts on the detection plate creates the interference pattern 
shown in Fig. [Tj 

The first calculation of the wave function double slit experiment using electrons^ was done 
using the Feynman path integral method.'^ However, this calculation has some limitations. It 
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FIG. 1: Schematic configuration of tlie experiment. 

covered only phenomena after the exit from the sUts, and did not consider reahstic shts. The 
sUts, which could be well represented by a function G{y) with G{y) = 1 for — /9 < y < P and 
G{y) = for \y\ > P, were modeled by a Gaussian function G{y) = e^^^/^^^. Interference was 
found, but the calculation could not account for diffraction at the edge of the slits. Another 
simulation with photons, with the same approximations, was done recently.'^ Recently, some 
interesting simulations of the experiments on single and double slit diffraction of neutron^ 
were done.'^ 

The simulations discussed here cover the entire experiment, beginning with a single source 
of atoms, and treat the slit realistically, also considering the initial dispersion of the velocity. 
We will use the Feynman path integral method to calculate the time-dependent wave func- 
tion. The calculation and the results of the simulation are presented in Sec. II. The evolution 
of the probability density of the wave packet just after the slits raises the question of the 
interpretation of the wave/particle dualism. For this reason, it is interesting to simulate the 
trajectories in the de Broglid^ and BohmP^ formalism, which give a natural explanation of 
particle impacts. These trajectories are discussed in Sec. III. 
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II. CALCULATION OF THE WAVE FUNCTION WITH FEYNMAN PATH IN- 
TEGRAL 



In the simulation we assume that the wave function of each source atom is Gaussian in 
X and y (the horizontal variables perpendicular and parallel to the slits) with a standard 
deviation cxq = = cxy = 10 fim. We also assume that the wave function is Gaussian in z 
(the vertical variable) with zero average and a standard deviation — 0.3 mm. The origin 
{x = 0,y = 0, z = 0) is at the center of the atomic source and the center of the Gaussian. 

The small amount of vertical atomic dispersion compared to typical vertical distances, 
~ 100 and 200 mm, allows us to make a few approximations. Each source atom has an 
initial velocity v = {vox,voy,voz) and wave vector k = {kox, koy, koz) defined as k = mw/h. 
We choose a wave number at random according to a Gaussian distribution with zero average 
and a standard deviation = = a^y = = may/\/?)h ~ 2 x 10^ m~^, corresponding 
to the horizontal and vertical dispersion of the atoms' velocity inside the cloud (trap). For 
each atom with initial wave vector k, the wave function at time t = is 

^o{x,y,z;kox,koy,koz) = ^o,{x;kox)^Oy{z;koy)^oAz'^hz) 

= (27r(T2)-i/4e-^V4-o%ifco.x 

X (27ra2)~i/^e-^'/^'^oeifco.y 
x{2-Kcxlr^/^e-''/^'''e""'^'. (1) 

The calculation of the solutions to the Schrodinger equation were done with the Feynman 
path integral method,'^ which defines an amplitude called the kernel. The kernel charac- 
terizes the trajectory of a particle starting from the point a = {xa,ya, Za) at time ta and 
arriving a at the point /? = {xp.yi^jZp) at time tp. The kernel is a sum of all possible 
trajectories between these two points and the times ta and tp. 

Using the classical form of the Lagrangian 

ip' ij^ z^ 

L{x,y,z,z,t)=m \- m \- m \-mgz. (2) 

2 2 2 

Feynman^^- defined the kernel by 

K{P,tf3;a,ta) ^ exp {-Sci{P,tf3;a,ta)) = exp {- L{x,y, z, z,t) dt), (3) 

with f^^ K{P, t/s; a, ta) dxa dya dza = 1. Hence, 

K{(3, tp; a, ta) = K^{xci, tp; a;„, ta)Ky{yp, tp] ya, ta)Kz{zi3, tp] Za, ta), (4) 

4 



with 



/ m im / {Xj3 — Xa^"^ 



2 



2 

y<exp — y^{zp + Za){t(3-t^) - ^{t(3-t^fj. (5c) 



For each atom with initial wave vector k, let us designate by to,; k) the wave function 
at time t^. We call S the set of points a where this wave function does not vanish. It is 
then possible to calculate the wave function at a later time at points (3 such that there 
exits a straight line connecting a and (3 for any point a G S. In this case, FeynmarP^I has 
shown that: 



iIj{(3, t; fcox, hy, hz) = / K{(3, t; a, to) ipia, ta] h^, hy, koz) dx^, dya, dza- (6) 

For the double slit experiment, two steps are then necessary for the calculation of the 
wave function: a first step before the slits and a second step after the slits. 

If we substitute Eqs. ([T]) and Q in Eq. ([6]), we see that Feynman's path integral allows 
a separation of variables, that is, 

^{x,y,z,t;kox,koy,ko^) = i)^{x,t; ko^)^y{y,t; koy)tp^{z,t; ko^). (7) 

References ITT] and treat the vertical variable z classically, which is shown in Appendix A 
to be a good approximation. Hence, we have z{t) = zo + vozt + gt'^/'2. The arrival time of 
the wave packet at the slits is ti(fo2, zq) = ^ ^(^i~^o) _[_ — For vqz = and zq = 0, 

we have ti = = 124 ms and the atoms have been accelerated to Vzi = gti = 1.22 m/s 
on average at the slit. Thus the de Broglie wavelength A = h/mVzi = 1.8 x 10~^m is two 
orders of magnitude smaller than the slit width, 2/im. 

Because the two slits are very long compared with their other dimensions, we will assume 
they are infinitely long, and there is no spatial constraint on y. Hence, we have for an initial 
fixed velocity voy-. 



'ipyiy, t; koy) = / Ky{y, t; = 0)ipo{ya; koy) dy^. (8) 
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Thus 

i^y{y.t;hy) = (27r3g(t))-^/^exp[- ^^ ~ ""^f) + ikoy{y - v^yt)], (9) 

with .oW=^o(l + ^). 

The wave packet is an infinite sum of wavepackets with fixed initial velocity. The proba- 
bility density as a function of y is 

r+oo 

py{y,t) = / {27rr')-'/'e-'^'y/'^'\^y{y,t;ky)\'dky = (2vreg(t))-V2 ,-.V2.gW^ (10) 



with el{t) = o"o(t) + (^)^ and a^lt) = o"o + (2^^)^- Because we know the dependence of the 
probability density on y, in what follows we consider only the wave function ipxixjt; ko^)- 



A. The wave function before the slits 



Before the slits, we have 



'ilJx{x,t;kocc) = (27rSo(t))~^/''exp [ + ikox{x - vqJ)] (11) 

Px{x, t) = (27re^(t))-V2 exp [ - ^] . (12) 

It is interesting that the scattering of the wave packet in x is caused by the dispersion of the 
initial position ctq and by the dispersion r of the initial velocity Vq^ (see Fig. [2]). Only 0.1% 
of the atoms will cross through one of the slits; the others will be stopped by the plate. 



B. The wave function after the slits 

The wave function after the slits with fixed zq and koz = mvoz/fi for t > ti{voz,zo) is 
deduced from the values of the wave function at slits A and B (cf.Fig. [s]) by using Eq. (|6]). 
We obtain: 

^px{x, t; fcox, koz, zo) =i)A + ^B, (13) 

with 



i^A= I Kx{x,t;Xa,ti{voz,Zo))ipx{xa,ti{voz,Zo);kox)dxa (14a) 

J A 

ipB = Kx{x,t;Xb,ti{voz,zo))ilJx{xb,ti{voz,zo);kox)dxb, (14b) 



B 
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FIG. 2: Evolution of the density px{x,t) before the shts, The time of the scattering of the initial 
wave packet {t is determined from z = zo + vozt + ^). 



where ipx{xa,ti{voz, Zq); ko^) and ipx{xb,ti{voz, Zq); ko^) are given by Eq. ([TT]) whereas 



Kx{x,t;Xa,ti{voz,zo)) and Kx{x,t;Xb,ti{voz,zo)) are given by Eq. ( [5a| . 
The probabihty density is 



Px{x,t;koz,Zo) 



-1/2 



exp 



- ^j\ilJx{x,t;kox,koz,zo)\'^dk 



Ox- 



(15) 



The arrival time t2 of the center of the wave packet on the detecting plate depends on zq and 
vqz- We have ^2 = 



a 



a ^ a 

and the atoms are accelerated to Vz2 = gt2 = 1.93 m/s. 



For Zq = and vqz = 0, t2 



2{h+l2) 



196 ms 



The calculation of Px{x, t; koz, zq) at any (x, t) with koz and zq given and t >ti is done by 



a double numerical integration: (a) Eq. (15) is integrated numerically using a discretization 



of kox into 20 values; (b) the integration of Eq. ( 13 ) using Eqs. ( 14 ) is done by a discretization 
of the slits A and B into 200 values each. Figure |4] shows the cross sections of the probability 
density (IV^a+V'-bP) for zq = 0, vqz = (/cqz = 0) and for several distances (A^ = ^gt"^ — ^gtl) 
after the double slit: 1pm, 10 p m, 0.1mm, 0.5 mm, 1mm and 113 mm. 




z 

FIG. 3: Schematic representation of the experiment: calculation method of the wave function after 
the slits. 

The calculation method enables us to compare the evolution of the probability density 
when both slits are simultaneously open (interference: {ipA + "^sP) with the sum of the 
evolutions of the probability density when the two slits are successively opened (sum of two 
diffraction phenomena: (iV'Ap + lV'Bp)- Figure [4] shows the probability density (|^/'aP + |^bP) 
for the same cases. Note that the difference between the two phenomena does not exist 
immediately at the exit of the two slits; differences appear only after some millimeters after 
the slits. 

Figures 5-7 show the evolution of the probability density. At 0.1mm after the slits, we 
know through which slit each atom has passed, and thus the interference phenomenon does 
not yet exist (see Figs. |4]and[7]). Only at 1mm after the slits do the interference fringes 
become visible, just as we would expect by the Fraunhoffer approximation (see Figs. |4] and 
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FIG. 4: Comparison between lipA + ips]'^ (plain line) and IV'Ap + IV'sP (dotted line) with ^0 = 
and koz = at (a) 1 /xm, (b) 10 m, (c) 0.1mm, (d) 0.5mm, (e) 1mm, and (f) 113mm after the 
slits. 
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FIG. 5: Evolution of the probability density Px{x, t; koz = 0,zq = 0) from the source to the detector 



screen. 



C. Comparison with the Shimizu experiment 



In the Shimizu experiment, atoms arrive at the detection screen between t = tmm and 
^max- To obtain the measured probability density in this time interval, we have to sum the 
probability density above the initial position zq and their initial velocity vqz compatible with 



^min ^ ^2 ^ ^maxj that is. 



Pxix, t2; koz, zo) e-^lJ-^-'e-^l-^-l dk^zdzo- (16) 

max 



The positions at the detection screen can only be measured to about 80 /im, and thus to 
compare our results with the measured results, we perform the average 

px+AO/im 



Pmeasured(-^) ^min ^ ^ ^ ^ri 



P{U, tmin < t < tmax) du. 



(17) 



80 fxm Jx-ioiim 

Figure |8] compares those calculations to the results found in Ref. [TTl The experimental 
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FIG. 6: Evolution of the probability density px{x,t; koz = 0,zq = 0) for the first millimeter after 
the slits. 



fringe separation is narrower than in our calculation, see Figs. |8]and Vi. This difference is 
explained by a technical problem in the Shimizu experiment.'^ 

III. IMPACTS ON SCREEN AND TRAJECTORIES 

In the Shimizu experiment the interference fringes are observed through the impacts 
of the neon atoms on a detection screen. It is interesting to simulate the neon atoms 
trajectories in the de Broglie-Bohm interpretation,'!^'!^ which accounts for atom impacts. 
In this formulation of quantum mechanics, the particle is represented not only by its wave 
function, but also by the position of its center of mass. The atoms have trajectories which 
are defined by the speed v(a;, y, z, t) of the center of mass, which at position (x, y, z) at time 
t is given bj» 



I ^S{x,y,Z,t) V log p{x,y,Z,t) X S h i ,*V7 l\ ^ T) ( l*V7 l\ ^ ^ 

\{x,y,z,t) = 1 = — [lm{ip vip) + Re(ip vip) x — 

m m mp |s| 



(18) 
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FIG. 7: Evolution of the probability density px{x,t; k{ 




slits. 



where ip{x,y,z,t) = a/ p{x, y, z, t) exp (^^(x, 2;, t)) and s is the spin of the particle. Let 
us see how this interpretation gives the same experimental results as the Copenhagen inter- 
pretation. 

If ip satisfies the Schrodinger equation, 



with the initial condition ^/^(a;, 2;, 0) = TpQ^x.y.z) = A/po(a^5 2/, -z) exp(|5'o(a;, 2;)), then p 
and S satisfy: 

with initial conditions S{x, y, z, 0) = So{x, y, z) and p(x, z, 0) = Pq{x, y, z). 

In both interpretations, p{x,y,z,t) = 2;, t)p is the probability density of the 

particles. But, in the Copenhagen interpretation, it is a postulate for each t (confirmed by 
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FIG. 8: Comparison of the probability density measured experimentally by ShimiziJ^ (left) and 
the probability density calculated numerically with our model (right). 

experience). In the de Broglie-Bohm interpretation, if pQ{x,y,z) is the probability density 
of presence of particles for t = only, then p{x,y,z,t) must be the probability density of 



the presence of particles without any postulate because Eq. (21) becomes the Madelung 
equation: 

f^+V- (pv) = 0, (22) 

(thanks to v = ^ + V x (^)s), which is obviously the fluid mechanics equation of conser- 
vation of the density. The two interpretations therefore yield statistically identical results. 
Moreover, the de Broglie-Bohm theory naturally explains the individual impacts. 

In the initial de Broglie-Bohm interpretation,'^^'^ which was not relativistic invariant. 



the velocity was not given by Eq. (18), but by v = — which does not involve the spin. In 
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the Shimizu experiment, the spin of each neon atom in the magnetic trap was constant and 
vertical: s = (0,0, |). In our case the spin-dependent term -^^^ x s = 2^(§^' 
neghgible after the sht, but not before. 

For the simulation, we choose at random (from a normal distribution /(O, 0, 0; 0"^, (Tk-, CTk) 
the wave vector k = {kox, koy, koz) to define the initial wave function ([T| of the atom prepared 
inside the magneto-optic trap. For the de Broglie-Bohm interpretation, we also choose 
at random the initial position (xo,?/0)^o) of the particle inside its wave packet (normal 
distribution /((O, 0, 0); (do, ctq, cTz)))- The trajectories are given by 

dx , , 1 dS h dp , 

— = v,{x,t) = -— + / 23a 

at m ox Imp oy 

dy , , 1 dS h dp / , X 

-£ = Vyix,t) = -— - — / 23b 
dt m oy Imp ox 

dz , . 1 dS , , 

— = Vz{x,t) = - — , 23c 
dt m oz 

where p{x,y,t; kox, koy) = \ipx{x,t; kox)'ipy{y,t; koy)\'^ and ip^ and ipy are given by Eqs. 
(fll. 



A. Trajectories before the slits 

Before the slits. Appendix B gives z{t) = zo^j^ + vozt + \gt^, x{t) = voxt + 
V^o + Vo^'^os (fit), and y{t) = Voyt + ^/xf+y^^ sin (p{t), with (p{t) = (fo + arctan( - 
TT^), COS090 = / ^° „ and sino^n = , ^° „ ■ For a given wave vector k and an initial 
position {xQ,yo, Zq) inside the wave packet, an atom of neon will arrive at a given position 
on the plate containing the slits. Notice that the term Vlogp x s/m adds to the trajectory 
defined by VS/m a rotation of — 7r/2 around the spin axis (the z axis). 

The source atoms do not all pass through the slits; most of them are stopped by the 
plate. Only atoms having a small horizontal velocity Vqx can go through the slits. Indeed an 
atom with an initial velocity vq^ and an initial position xq, yo arrives at the slits at t = ti at 
the horizontal position x(ti) = Vo^ti + yoicroiti) / ao) . For this atom to go through one of the 
slits, it is necessary that |a;(ti)| < x, with x = {d + b)/2 = 4 x 10"^ m and —2ao < yo < 2(Tq 
and ti ~ 0.124 s. Consequently, it is necessary that the initial velocities of the atoms satisfy 
\vqx\ < (a^ + '^(yo{ti)) /ti = Vqx — 3.9 X 10"^m/s. The double slit filters the initial horizontal 
velocities and transforms the source atoms after the slits into a quasi-monochromatic source. 
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slit A slit B 



slit A slit B 



FIG. 9: Trajectories of atoms before the slit. Note that if the initial velocities |fox| > 3.9x10 '^m/s, 
then no atoms can cross the slit. 

The horizontal velocity of an atom leads to a horizontal shift of the atom's impacts on the 
detection screen. The maximum shift is Ax = Vq^AI, where At is the time for the atom to 
go from the slits to the screen (At = ^2 — — 0.072 s); hence Ax ~ 2.8 x 10^^ m. This shift 
does not produce a blurring of the interference fringes because the interference fringes are 
separated from one another by 25 x 10^^ m ^ Ax. Note that if the source was nearer to 
the double slit (for example if li = 5 mm, then Ax ~ 10 x 10^^ m), the slit would not filter 
enough horizontal velocities and consequently the interference fringes would not be visible. 

The system appears fully deterministic. If we know the position and the velocity of an 
atom inside the source, then we know if it can go through the slit or not. Figure [9] shows 
some trajectories of the source atoms as a function of their initial velocities. Only atoms 
with a velocity \vox\ < vox can go through the slits. 

B. Velocities and trajectories after the slits 

In what follows, we consider only atoms that have gone through one of the slits. After 
the slits, we still have z{t) = vozt+ \gt^ + zq^cj z{t) / a z) , but now Vxit) and Vy{t) and x(t) and 
y{t) have to be calculated numerically. The calculation of Vx{x,t) is done by a numerical 
computation of an integral in x above the slits A and B (see Appendix B); x(t) is calculated 
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with a Runge-Kutta method.'^ We use a time step At which is inversely proportional to the 
acceleration. At the exit of the slit, At is very small: At ^ 10^® s; it increases to At ~ 10~^ s 



at the detection screen. Figure [10] shows the trajectories of the atoms just after the slits; Xq 
and yo are drawn at random, = 0, with vqx = "^oz = 0. 




C. Impacts on the screen 

We observe the impact of each particle on the detection screen as shown by the last 



image in Fig. 11 The classic explanation of these individual impacts on the screen is the 
reduction of the wave packet. An alternative interpretation is that the impacts are due to 
the decoherence caused by the interaction with the measurement apparatus. 
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FIG. 11: Atomic impacts on the screen of detection. 



In the de Broglie-Bohm formulation of quantum mechanics, the impact on the screen is 



the position of the center of mass of the particle, just as in classical mechanics. Figure 11 



shows our results for 100, 1000, and 5000 atoms whose initial position (xo,?/0)^o) are drawn 
at random. The last image corresponds to 6000 impacts of the Shimizu experiment. -"^^^ The 
simulations show that it is possible to interpret the phenomena of interference fringes as a 
statistical consequence of a particle trajectories. 



IV. SUMMARY 



We have discussed a simulation of the double slit experiment from the source of emission, 
passing through a realistic double slit, and its arrival at the detector. This simulation is 
based on the solution of Schrodinger's equation using the Feynman path integral method. A 
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simulation with the parameters of the 1992 Shimizu experiment produces results consistent 
with their observations. Moreover the simulation provides a detailed description of the 
phenomenon in the space just after the slits, and shows that interference begins only after 
0.5 mm. We also show that it is possible to simulate the trajectories of particles by using 
the de Broglie-Bohm interpretation of quantum mechanics. 



APPENDIX A: CALCULATION OF ipz{z,t;h 



Oz) 



Because there is no constraints on the vertical variable z, we find using Eqs. (g and Q 
for alH > (before and after the double slit) that 

ilj^{z,t;koz) = / Kz{z,t;zo,,t = 0) xi!o^{zo,;kQ-,)dZa, (Al) 
Js 

where the integration is done over the set S of the points Za where the initial wave packet 
■ipo^^Za, koz) does not vanish. We obtain: 



im 



ij,{z,t;koz) = (27rs2(t))-i/4gxp 
X exp 

where Sz{t) = az{l+ 2Wj)- Consequently we have 
witha.(t) = |s.(t)| = a4l+(^)y/l 



{{voz + 9t){z-Vozt/2)-'^^) 



(A2) 



[z-vo,t-gt y2)^- 
2al{t) 



Note that <Jz(t) is negligible compared to li (a^ = 0.3mm and 



ht 



(A3) 



10 ^ mm are 



negligible compared to h = 76 mm for an average crossing time inside the interferometer of 



t ~ 200 ms). Therefore {2'Ka^^{t))-^^^ exp [ 



{z-vo^t-gt^/2f 



] ^ 6o{z - vozt - 3^), and if zo 



IS 



the initial position of the particle, we have 2; ~ + vozt + ^ at time t. 



APPENDIX B: CALCULATION OF THE ATOM'S TRAJECTORIES 



The velocity (18) applied to Eq. (A2) gives the differential equation for the vertical 
variable z, 

dz , , 1 dS {z-VQzt-gt^l2)hH 

- = Vz{z, t) = -— = v,z + gt+ ^ , 2 2 ' ^1 

dt m oz Am^alal[t) 
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from which we find 



z{t) = Vo,t+-gt^ + Zo^^. 

Z cr. 



(B2) 



Equation (B2) gives the classical trajectory if = (the center of the wave packet). 

The velocity ( 18 ) applied before the slit to Eqs. ^ and ( 11 ) gives the differential equations 
in the x and y directions: 

^ = v^{x,t) = 

at m ox Imp oy 

dy 1 dS h dp 



dt 



1 OS h dp (x — voxt)hH h{y — voyt) 

+ — = Vqx + 



m dy 2mp dx 



Voy 



{y-Voyt)hH h{x-Vo,j) 



+ 



(B3a) 
(B3b) 



It then follows that 



xit) = voxt + + I/q-^ cos (fit) 

V (To 

y{t) = voyt+Jx^ + y^ smifit), 

V (To 



(B4a) 
(B4b) 



withv9(t) = v5o + arctan(-^^), cos(v?o) 



= , and sin(y9o 



yo 



f. Equations (B4a) 



and (B4b) give the classical trajectory if xq = yo = (the center of the wave packet). 



After the slits, the velocity Vx{x^t) — ^ 
using Eqs. (|6|, (14a) and (14b). We obtain: 



given by Eq. (18) can be calculated 



Vx{x,t) 



t-ti 



X + 



2(a= + Pf) 



H{x,t) 



H{x,t) 



with 



H{x,t) = / f{x,u,t)du+ / f{x,u,t)du 

JxA-b JxB~b 

cix,t) = [fix,u,t)]:zifj,+[f{x,u,t)]:ziiti 

where Xa and Xb are the centers of the two slits, and where 

f{x, u, t) = exp [{a + iPtW + hx,M 
1 

a = 



Pt 



7x,t 



m 
2h 



t - h 



mx 



h 1 



Ml 



h{t - ti 



- Ptlx,. 
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APPENDIX C: CONVERGENCE 



In simulation, it is possible to make the Planck constant h tend towards zero. So, it is 
possible to show, about the example of the Young slits, the convergence of the quantum 
mechanics to the classical mechanics. In the following simulations, one take kox = koz = 0. 




FIG. 12: Cross section of the probabiUty density on the screen of detection when h is divided by 
5, 10 and 100. 



The figure 12 shows, when h is divided by 100, how the interference fringes are getting 
strongly more and more narrow up to the distance between the slits. 



The figure 13 shows, when h is divided by 1000, that the interference fringes disappear 
and the probability density converges on the classical probability density {h = 0). 



The figure 14 shows how the trajectories become strongly narrower when h is divided by 



100. 
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FIG. 13: Cross section of the probability density on the screen of detection when h tends towards 



zero. 



The figure 15 clearly shows the convergence of the quantum trajectories to the classical 
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FIG. 14: Trajectories of atoms for whole experience when h is divided by 5, 10 and 100 



trajectories (/i = 0). 

After showing results of the evolution of the probability density and of the trajectories 
of atoms, it is difficult to throw down the hypothesis of trajectory. 

Everything seems to happen also as if atoms have well a trajectory. 
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FIG. 15: Trajectories of atoms for whole experience when h tends towards zero. 
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